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i We develop a general formalism to investigate the effect of quenched fixed charge disorder on 

effective electrostatic interactions between charged surfaces in a one-component (counterion-only) 
Coulomb fluid. Analytical results are explicitly derived for two asymptotic and complementary 
cases: i) mean-field or Poisson-Boltzmann limit (including Gaussian- fluctuations correction), which 
^v^j . is valid for small electrostatic coupling, and ii) strong-coupling limit, where electrostatic correlations 

■ mediated by counterions become significantly large as, for instance, realized in systems with high- 

valency counterions. In the particular case of two apposed and ideally polarizable planar surfaces 
with equal mean surface charge, we find that the effect of the disorder is nil on the mean-field level 
and thus the plates repel. In the strong-coupling limit, however, the effect of charge disorder turns 
out to be additive in the free energy and leads to an enhanced long-range attraction between the 
two surfaces. We show that the equilibrium inter-plate distance between the surfaces decreases 
for elevated disorder strength (i.e. for increasing mean-square deviation around the mean surface 
charge), and eventually tends to zero, suggesting a disorder- driven collapse transition. 
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O ; I. INTRODUCTION 

o . 

Electrostatic interactions usually provide the repulsion that stabilizies charged colloids and are one of the two 
essential ingredients in the DLVO theory of colloidal stability 0, 0. Electrostatic interactions in the presence of 
pi , mobile counterions are standardly described by the Poisson-Boltzmann theory Q, Q that embodies the mean-field 
pq ■ approach to Coulomb fluids and leads to pronounced repulsive interactions between like-charged macroions. It has 
recently been realized however 0, IE IE S H 13, HO, HH Hi, EE HI, EE EE HI, HE IH, E3, EH [23, HE IM [2E HE p3l, 
that in the presence of polyvalent counterions, electrostatic interactions can mediate strong attractive interactions 
between like-charged macroions. This attraction can not be captured by the mean-field approach and a new paradigm 
dubbed the strong-coupling limit |2 ll |25| was devised to describe the equilibrium properties of Coulomb fluids when 
the mobile counterion charges become large. The transition from the mean-field Poisson-Boltzmann description to 
the strong-coupling limit is governed by a single dimensionless electrostatic coupling parameter, being a ratio of two 
lengths, namely, the Bjerrum length, which identifies Coulombic interaction between ions themselves, and the Gouy- 
Chapman length, which describes electrostatic interaction between the ions and the charged macroion surface. This 
i-^J ■ ratio involves the char ge v alency of counterions, the dielectric constant of medium and the surface charge density of 
interacting macroions |27| . The emerging picture of equilibrium properties of Coulomb fluids has thus become much 
O ■ richer than conveyed for many years by the DLVO paradigm. 

The strong-coupling attraction is not just a refinement of the mean-field description. In fact it reverses the role 
J> , of Coulomb interactions in the DLVO theory [27|. Instead of stabilizing the charged macroions they act more like 
Lifshitz-van der Waals interactions that tend to collapse them. Since Coulomb interactions in the strong-coupling 
limit are much stronger than the Lifshitz-van der Waals interactions, they themselves govern the destabilizaition of the 
colloids. The collapse of a highly charged polyelectrolyte, such as DNA, in the presence of polyvalent counterions is 
the most drammatic example of unexpected and counter- intuitive features of the strong-coupling electrostatics 1 11.1251. 
Though other physical mechanisms, such as hydration layer complementarity certainly contribute fine details to 
this collapse and account for its molecular specificity, electrostatic correlations involved in the strong-coupling limit 
provide the universal background for this intriguing phenomenon. The elucidation of this collapse in terms of the 
strong-coupling electrostatic interactions represents one of the major recent advances in the theory of colloid stability. 
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In what follows, we shall concentrate on yet another facet of Coulomb interactions in charged colloids, namely, the 
quenched disordered distribution of surface charges. Usually the charge on the surfaces of the macroions on the Poisson- 
Boltzmann or the strong-coupling level is assumed to be homogeneous and constant. This is in general quite a severe 
assumption and there are well known cases where this assumption is not realistic at all |28| . Random polyelectrolytes 
and polyampholytes present one such case |2j|, |3(J, |3l| . There the sequence of charges can be distributed along the 
polymer backbone in a disordered manner. The Coulomb (self-)interactions of such polyelectrolytes are distinct and 
different from homogeneously charged polymers. 

A case even closer to the present line of reasoning are investigations of interactions between solid surfaces in the 
presence of charged surfactants. The aggregation of surfactants at solid surfaces in aqueous solutions was investigated 
with atomic force microscopy and shows structures consistent with half-cylinders on crystalline hydrophobic substrates 
for quaternary ammonium surfactants (above the critical micelle concentration), full cylinders on mica, and spheres on 
amorphous silica [32j,|33|. Such interfacial aggregates whose emergence and structural details depend on the method of 
preparation result from a compromise between the natural free curvature as defined by intermolecular interactions and 
the constraints imposed by specific surfactant-surface interactions and can pattern interacting surfaces at nanometer- 
length scales. Similar interfacial structures are also seen for interacting hydrophilic mica surfaces in the presence 
of cetyl-trimethyl-amonium bromide (CTAB) 34]. The surfaces appear to be covered by a mosaic of positively and 
negatively charged regions and experience a strong, long-ranged attraction which is comparable in magnitude to that 
between hydrophobic surfaces, and is orders of magnitude larger than that expected from Lifshitz-van der Waals 
forces [34|. The patterning of interacting surfaces described above is highly disordered, depends on the method of 
preparation and has basic implications also for the forces that act between other types of hydrophilic surfaces with 
mixed charges, most notably the surfaces of cells and proteins, as well as in synthetic systems. 

Motivated by these observations we set ourselves to investigate the effect of quenched disordered surface charge, 
specified by a constant mean value and mean-square deviation, on the interactions between charged surfaces in ionic 
solutions. We formulate the general partition function of a system of charged objects with a Gaussian-distributed 
quenched surface charge in a one-component Coulomb fluid using field-theoretical techniques and the standard replica 
trick. As a particular case, we focus on the effective interaction between two apposed, ideally polarizable planar 
surfaces of equal mean and mean-square deviation for disordered surface charge density. An explicit form for the 
interaction free energy can be obtained on the Poisson-Boltzmann level as well as on the strong-coupling level. The 
former, interestingly enough, shows no effect of the quenched charge disorder. The interaction free energy at the 
mean-field level plus quadratic (Gaussian) fluctuations around the mean-field solution is exactly the same as the one 
for a uniform charge distribution characterized solely by a constant mean surface charge. In the strong-coupling limit, 
the situation is altogether different. Here the quenched surface charge disorder contributes an additive attractive tail to 
the interaction determined by the mean surface charge. This attractive contribution is linear in the disorder coupling 
strength, which is proportional to the mean-square deviation of the charge distribution, and is a direct analogon of the 
electrostatic coupling parameter introduced in the case of no charge disorder |27| ■ This attractive tail can completely 
overwhelm the entropic repulsion due to counterion confinement on the strong-coupling level and can thus completely 
eliminate the finite equilibrium distance between the surfaces in this limit |24|. In other words, the introduction of 
charge disorder can destabilize the closely-packed bound state of two uniformly and like-charged surfaces [24l [27j , 
suggesting a collapse transition. This is an important result because especially in the biological context the surface 
charge disorder is ubiquitous, stemming not only or not at all from the method of preparation of the macroions, but 
from their intrinsic molecular disorder, dictated eventually by the evolutionary processes. 

The outline of the paper is as follows: First we enunciate the results, which is indeed quite simple, but can only be 
derived by fairly detailed and rather involved technical arguments. The general formulation of the partition function 
for a one-component Coulomb fluid between charged planar surfaces is developed next in the form, where the trace 
over the quenched disorder has already been evaluated exactly on the replica level. This very complicated partition 
function in the replica space can then be explicitly traced over local electrostatic field fluctuations on the Poisson- 
Boltzmann level as well as on the strong-coupling level. This can be accomplished by noting the properties of a 
certain type of symmetric matrices. The formalism developed in the present work is quite general and can be applied 
to investigate charge disorder effects in a variety of systems, including interactions between spherical and cylindrical 
macroions, as well. 

II. MAIN RESULTS 

Since the derivation of our main result contains several rather subtle and highly non-trivial technical points we 
enunciate it first stripped of all technicalities. Our model system is composed of two ideally polarizable surfaces 
located at z — ±a, with the inter-surface spacing D — 2a (see Fig. The mean surface charge density on both 
surfaces is the same with a magnitude equal to cr, while the mean-square deviation in surface charge density due 
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FIG. 1: Schematic representation of our model system: two ideally polarizable planar surfaces with quenched surface charge 
disorder at separation D with mobile (point-like) counterions of charge Zeo confined in the inter-surface space. The surface 
charge density is characterized by a Gaussian distribution with fixed mean value of magnitude a and mean-square deviation g. 



to quenched disorder is given by g. The one-component Coulomb fluid between the two surfaces is composed of 
counterions of charge Ze , where Z is the counterion valency and eo is the elementary charge. For uniformly charged 
surfaces (g — 0), the electrostatic coupling parameter that identifies the strength of electrostatic correlations in the 
system is given by (27J 



S = 2nZ 3 £ 2 B a = (1) 
A* 

where £b = ^^/(ATTeeo) is the Bjerrum length (with f3 — 1/fcrjT) and \i = eo/ (2it£bZ<t) is the Gouy-Chapman length 
associated with the mean charge. As was demonstrated by Netz |24| for the undis ordered case, the partition function 
of the one-component Coulomb fluid has two different fixed points characterized by the value of S. For S <C 1, the 
partition function is given by the Poisson-Boltzmann (mean-field) fixed point. Fluctuations around this fixed point 
may be accounted for by a systematic loop expansion 0, H, El 0] • At elevated coupling, the loop expansion scheme 
breaks down and the partition function for S ^> 1 is determined by the asymptotic strong-coupling theory, which 
corresponds to the leading term of a virial expansion, and is in fact given by single-particle contributions. 

In the presence of quenched disorder (<? > 0), one can follow the same line of reasoning, by introducing the mean 
electrostatic coupling parameter 3 in analogy with Eq. Q but with a now being the mean surface charge density. 
In the mean-field limit, characterized by 5 < 1, we demonstrate that the disorder makes no contribution to the 
interaction free energy, if the bounding surfaces are ideally polarizable and thus do not allow the penetration of the 
field. The total interaction pressure in this limit is given by the standard mean-field (Poisson-Boltzmann) expression 
Pmf(D) and an additive contribution, P2(D), due to quadratic (Gaussian) fluctuations around the mean-field solution 
as 

P{D)=P MF {D)+P 2 {D). (2) 



Both terms depend only on the mean surface charge density a and have been thoroughly analyzed before |7], llCt lla 
liH I24I ] . The Poisson-Boltzmann limit for ideally polarizable bounding interfaces is thus unaffected by the presence of 
disorder, i.e. it does not depend on g. This result is completely consistent with conclusions reached in Ref. |28j| via 
an altogether different route. 

In the strong-coupling limit, characterized by large mean electrostatic coupling parameter H » 1, however, the 
interaction free energy is obtained to contain an essential contribution stemming from the quenched surface charge 
disorder. The total interaction pressure is given by the sum of the standard strong-coupling pressure Psc(D) [24j 
that depends only on the mean charge density a plus a disorder contribution, P X (D), depending on g, i.e. 

P(D)=P SC (D)+P X (D), (3) 

where the disorder term in dimensionless representation reads 

m =-<!> 
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which is attractive and linear in \ defined as 

Z 2 (f3e ) 2 g 22 

x = "W = 2nZ (5) 

Clearly by comparing Eq. (yi and Eq. J3J), we can claim that x represents the dimensionless disorder coupling 
parameter. Psc(-D) has been analyzed thoroughly before [24^. The disorder term in Eq. (J3J) incidentaly has the 
same scaling form as the one-particle ideal entropy term in Psc(D). It can change the sign of the interactions at 
small values of the inter-surface spacing £>, as will be shown later. The quenched charge disorder thus affects only 
the strong-coupling limit in a rather simple way, contributing an additive term to the free energy or the interaction 
pressure. After this introductory exposition, we next derive these results in full glory and all the relevant technical 
details. 

III. GENERAL FORMALISM: THE REPLICA METHOD FOR QUENCHED AVERAGING 

The partition function of a one-component (counterion-only) Coulomb fluid in the field of an external fixed charge 
charaterized by the volume charge density p(r) can be derived in the form of a functional integral over the fluctuating 
electrostatic field 0(r) as 0, HI 

z = e -ilndct/Mr,r') f £)[0( r )] e -^H[*(r)] ) (g) 

where the field-action is given by 

pn[<j}(r)} = y I /drdr' ^(r)^ 1 (r, r')^(r') - A fdr fi(r) e' lf3Ze ° <l ' {r) + i(3 [dr p(r)0(r). (7) 



Here the geometry function f2(r) specifies the volume accessible for counterions in space (it, for instance, takes into 
account the presence of hard walls). The rescaled absolute activity (fugacity), A (connected to the chemical potential, 
/i, via fj, = In A), is defined as 

A = Ae5 e o/M r ' r ). ( 8 ) 
The Coulomb interaction potential v(r,r') = l/(47reeo|r — r'|) is a solution of 

-eeoV 2 f(r,r') = <5(r-r'). (9) 
and thus obviously the inverse Coulomb operator is given by 

u -1 (r,r') = -ee V 2 £(r - r'). (10) 

We assume in what follows that the external charge distribution has a quenched disordered component and is in fact 
distributed with a Gaussian probability distribution around its mean po(r) as 



const, x e 2 



Here the width of the charge disorder distribution, i.e. the mean-square charge density is given by q(r). The average 



:Y Tl 
O ii 



over quenched charge disorder is obtained by applying the standard Edwards- Anderson replica ansatz [35j in the form 

Z"- 1 



T = -k B T logZ = -k B T lim , (12) 

n^o n 

where the disorder average is defined with respect to the external charge density distribution, Eq. i|ll[l . as 

(T7)= /p[p(r)](...)e-5/ dr ^ 1 (r)(p(r)-p (r)) 2 ^ (13) 
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The quenched average over the charge density disorder affects only the source term in Eq. Q . Thus we only need to 
evaluate 

e - l /3/dr P {r)4>(r) = J V[p(r)} e ~ i(3 ^I dr P( r ) <f>a 0) - 3 / dr g~ 1 (r) O(r)-po) 2 = 

= const, x e~^ 2 / dr - 9(r) M')M r ) ~ # E« I dr Po(') E« <M0 

(14) 

for a, P — 1, . . . ,n being the replica labels. Taking this into account, the final form of the replicated partition function 
Z n can be obtained as 

-^ =e -i„lndet/Mr,r') J V ^ a ( r )] e ~i3H[4> a (p)] > ( 15 ) 

with 

/8W[^a(r)] = 11]// drdr ' ^M*, r')0a(r)^(r) - 7[0„(r)] + i/3 / dr Po (r) £ a (r), (16) 

a/3 " a 

where 

Vfo a (r)] = A /dr fi(r) ^e^ Ze ^»W. (17) 

a 

Above we introduced the following matrix in the replica space 

V a0 (r, r') = (/JiT^r, r') + (3 2 g(r)S(r - r')) <W + (3 2 g(r)S(r - r')(l - <M- (18) 

The expression (11511 together with Eq. 112(1 represents the starting formulation for the free energy in the presence of 
quenched charge disorder. There are two separate problems with the above replicated partition function l|15|) that we 
have to resolve. First of all there is the non-linear term in the field-action that precludes direct integration. On top 
of that we have to evaluate the functional integral in the replica space for arbitrary number of replicas, n, and then 
take the limit of n — ► 0. 

The way we will approach this rather formidable task is to combine the methods developed for the one component 
(counterion-only) Coulomb fluid without disorder [24J and modify them as we proceed to incorporate appropriately 
the disorder effects. We shall start with the Poisson-Boltzmann or the saddle-point limit, which is exact for S — > 
p4|. and then proceed to the strong-coupling limit. 



IV. THE MEAN-FIELD (POISSON-BOLTZMANN) LIMIT 

Let us first investigate the saddle-point limit of the functional integral (|15fl . The saddle-point equation for the 
i> a (r) field is given obviously by 



/ dr'X> Q/3 (r,r')<Mr') - +i{Jf>o(r) = 0. (19) 

Taking into account Eq. (|18|l . this amounts to 



-/? eeo v 2 0a(r) + P 2 g{v) M*) - dV }^ ] + = °- ( 2 °) 

p (#a(r) 

In the replica formulation we have to take the limit n — * 0, which implies that 

(21) 

n— »0 * — 4 



By furthermore making the substitution <fi a (r) — > i<p a ( r ) an d noting that in the limit n — > 0, the index a becomes 
irrelevant, and denoting the corresponding potential as 0mf(f), we are left with 



ee V 2 ^»MF(r) + \Ze fi(r) e^W') - -p (v). (22) 
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This is of course nothing but the Poisson-Boltzmann equation for the fixed mean charge density po( r ) and the mean- 
field Poisson-Boltzmann potential </>MF( r ) lEHS- The corresponding mean- field (saddle-point) action can be written 
as 

(3n[cf>M F (r)} = -ieeo/3 J dr (V0 MF (r)) 2 - A J dr fi(r) e-"^W + /? / dr jOo (r)0 M F(r). (23) 

The effect of the quenched charge disorder on the mean potential is thus nil. Obviously the saddle-point action can 
be taken out of the functional integral since it is simply multiplicative in the number of replicas. 

We now discuss the fluctuations around the mean-field, or the saddle-point, fields. Expanding f3H{(j) a (r)} in Eq. 
(|f 6[) to the second order in the deviations from </>mf (r) , we remain with 



where 7t£mf [0mf (r)] = W[0 M F(r)], 



/37i 2 [0 Q (r)] drdr'0 o (r)a ofl (r,r')^(r), 

Q/3 



and 



£ Q/3 (r,r') = £> Q/3 (r,r') 



a 2 y[^ Q (r)] 



<5 Q( 3<5(r - r'). 



9^(r) 

In the above expression the second term has to be evaluated at the saddle-point and thus yields 

3 2 V[«lo(r)] 



-A(/3Ze ) 2 ft(r) e -0S»°*M*M. 



(24) 



(25) 



(26) 



(27) 



0MF 



Thus we obtain explicitly 

g a p{v,v')= [^- 1 (r,r') + (A(^e ) 2 fl(r) e^^W + /3 2 5 (r)) <5(r - r')] <W + /3 2 g(r)5(r - r')(l - M- (28) 

Since ,£W2[0a(r)] is Gaussian in <f> a (r) the functional integral, Eq. 112 111 , can be evaluated explicitly yielding 

-gn_ e -inlndet^u(r, r') - n/3WMF[<frvtF(r)] - 5 lndet £ Q| o(r, r') , 2g s 

The task now is to evaluate the determinant det £ Q/ g(r, r') explicitly as a function of n and then use this to evaluate 
the n — > limit. This can be done with the help of the following matrix identity. Take a symmetric n x n matrix of 
the general form 



M al3 = 



b s s 
s b s 
s s b 



= b S af 3 + s (1 - 6 a p). 



(30) 



It is rather straightforward to show via induction that 



det M al3 = (b- s) n 1 



b- s 



In order to exploit this matrix identity we first define the operator G(r,r') as 

dr" Uv- 1 ^,*") + \(f3Ze ) 2 fi(r) e-^^W^r - r")) G(r",r') = 5(r-r'), 



or equivalently 



(-/3ee V 2 + \{f3Ze ) 2 fi(r) e" 



/3Ze o 0MF(r) 



G(r,r') = 5(r-r'). 



(31) 



(32) 



(33) 
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With this definition and by using the Trlog formula, the determinant det G a p{r, r') can be obtained explicitly as 

i In det g aff (r, r') = -f In det G(r, r') + ±Tr In (S(r - r') + n/3 2 5 (r)G(r, r')). (34) 
The final expression for the replicated partition function becomes 

lnS" = -|hidet Jdx' Pv(r,r')G- l (r' ,r") - n/37Y M F^MF(r)] - ±Tr In (<S(r - r') + n/3 2 5 (r)G(r, r')) . (35) 

We have thus evaluated the partition function on the Poisson-Boltzmann level plus Gaussian-field fluctuations exactly 
as a function of n. The replica trick now leads directly to the free energy, Eq. (|12fl . of the form 

Z" - 1 

T\ = —k B T lim = 

n— >0 n 

= Wmf^mfW] + ^Indct /dr'/MryjG-V,!-") + §Tr 3 (r)G(r,r'). (36) 

These three terms have a straightforward physical interpretation 0, . The first term is the standard mean-field 
result, the second is the one-loop correction around the mean-field, and the third one is the contribution of the 
disorder. Obviously the latter is linear in the strength of the disorder g. 

We now go back from the grand-canonical to the canonical partition function via the Legendre transform p4| 

Fn =Tx + k B TN In A (37) 

where N is the total number of counterions, which can be expressed through the absolute activity, A, as 

N = -X d -^ = -X d -m. (38) 

(This simply means that the log of the absolute activity is defined up to an additive constant.) To the lowest order, 
and on the saddle-point level, we obtain straightforwardly 

N = A / dr 0(r) r" a ^ MF(r) (39) 



This relation is nothing but the normalization condition for the mean-field density profile of counterions, 
A fi(r) e - ?Ze »* OT W, that can be connected with the electroneutrality condition via 2a S = NZeo, where S is the 
total area and a the surface charge density of macroions. Thus we finally remain with 

T N =F\- k B T Nlnjdr fi(r) e -^ e »^W + k B T N\nN. (40) 

For two charged planar surfaces, for instance, the interaction pressure can be obtained by differentiating the free 
energy with respect to the separation, D. 



V. MEAN-FIELD RESULTS FOR TWO CHARGED PLATES 



Let us assume now that our system is composed of two planar surfaces, located at z = ±a, with mean surface 
charge density — a. It is thus obviously homogeneous in directions p — (x,y). The mean volume charge distribution 
therefore reads 

po(r) = po(z, p) = -aS(z — a) - <rS(z + a). (41) 
The mean-square charge density due to disorder is analogously written as 

ff(r) = g(z, p) = gS(z-a)+ gS(z + a). (42) 
Note also that since counterions are confined in the inter-surface space (Fig. QJ, the geometry function fi (r) reads 

!1 —a < z < a 
(43) 
otherwise. 
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From Eqs. i|36|) and (|40|l we derive the corresponding canonical free energy in the form 

Fn =T N {g = {)) + §g[d 2 p (G(a, a; p, p) + G(-a, -a; p, p)) , (44) 



where we decoupled the disorder independent terms, ^Nig = 0), from those that depend on g. The former contribution 
may be written as J-N(g = 0) = .F^P^mfM] + fffi, where 

•^ F [</>MF(r)] = /dr (-^(V0 MF (r)) 2 +po(r)<Wr)) ~k B T Nln J dr Q(r) e-^ Ze ^^+k B T (NlnN-N). (45) 

is the mean-field Poisson-Boltzmann free energy, and J-ffi = In det J dr'/3w(r, r')G _1 (r', r") is the contribution of 
Gaussian fluctuations around the mean-field solution. These terms have already been evaluated elsewhere Q and will 
not be analyzed here. 

If we assume furthermore that the mean electrostatic field is confined solely to the region between the surfaces, 
that is assuming that the surfaces are ideally polarizable, then the last term in Eq. (|44J) describing the effects of the 
disorder is identically zero since it is straightforward to show that in this case 

G(a, a; p, p) = G{— a, —a; p, p) = 0. 

Therefore, with ideally polarizable bounding surfaces, the effect of the disorder on the electrostatic interactions is zero 
also on the one-loop level. Evaluating the interaction pressure between the two apposed surfaces from Eq. (|4^|) . we 
thus find 

P(D) = Pmf(D) + P 2 (D). (46) 

The first component of pressure is obtained from the mean-field Poisson-Boltzmann free energy Ij45(l , which is always 
repulsive. For two planar equally charged surfaces at large separations and in dimensionless units, one has 0,0] 

HD) = P^*(%Y. (47) 



2nl B a 2 V D 

The second term, stemming from fluctuations, in fact corresponds to zero- frequency (static) Lifshitz-van der Waals 
interaction 0- For two plates, this term contributes an attractive component that scales as ylEElElEI 

tW-&%*^(£\'uPIA- (48) 



2n£ B a 2 \D 

The preceding result is pleasingly consistent with a previous study of quenched disorder on the weak-coupling level 
|28j |. which also shows no disorder contribution to the density profile of counterions at a single charged surface. 

VI. THE STRONG-COUPLING LIMIT 

Now we address the strong-coupling limit, S > 1, which was shown to correspond to the first order term in the 
expansion of the partition function in terms of the absolute activity |2d| . Expanding Eq. (|15|l in terms of A we obtain 
to the lowest order 

= 2% + A Zf + 0{\ 2 ) = 



e -i n lndet/Mr,r') p[0 Q ( r )] e -/3Wo[^W] + \ Q( R ) J^J V [^)\ e^^-M'*^ 



(49) 



with 



/^o[i(r)] = iW/ drdr'^(r,r')0 a (r)0 /3 (r)+0 J dr p (r) ]T <f> a (r), (50) 

a/3 a 
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and 

0Wi[0 a (r)] = ± Wf drdr' V*fl{ry)4> a (r)M r ) -*0 I dr Po(r) 5> Q (r) + i/?Ze o 7 (R). (51) 

T> a p(r, r') was defined already in Eq. JTSJ|. The first functional integral in Eq. (|49|l can be evaluated straightforwardly 
by noting that it is equivalent to Eq. (|15|l for V[0 a ] = 0. Since in this limit obviously 

(3G(v,v')=v(r,r') 

we obtain immediately 

7, n — 1 



J- a = — fceT lim 

n— >o n 

= i ^ drdr' Po (r)«(r, r>o(r') + iTr g(r)v(r, r') (52) 

corresponding to the zeroth-order term of the virial expansion. The first term is the standard electrostatic interaction 
energy between external charges and the second term is the contribution of the disorder. We note here again that 
with the ansatz ()42J) . the disorder term does not depend on the separation between the bounding surfaces because 
v(z = a, z' = a : p, p) does not depend on a. 

The second term in Eq. I|49[l is a bit trickier but can be again straightforwardly reduced to 

/n 
dRfi(R)^e^ (R) , (53) 



7=1 



where 



(3u y (R) = (3 2 Ze Q f dr'^I^r'.RVoCr') - 7}(f3Ze ) 2 T>~^(R, R). (54) 

The task now is to evaluate the expressions ^2 a P^(r', R) and £>~^(R, R) explicitly as a function of n. This can be 
done just as before by considering the following identity valid for a symmetric n x n matrix of the general form 

M al3 = bS a0 + s (1-S a0 ). (55) 

One can derive via induction that 



and 





With the help of these matrix identities we can derive that 



(R, R) = J dr G(R, r)T(r, R), (58) 

where 

T(r, r') = «5(r - r') - (3 2 J dr" 5 (r")G(r", r) (S(r - r') + n(3 2 g (r)G(r, r')) _1 . (59) 

Similarly 

^I?- 7 1 (r',R) = |dr'G(r,r')(<5(r-r')+n/3 2 5(r)G(r,r')) _1 . (60) 
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Neither of the above two equalities depends on the index 7, and the 7-terms in Eq. <|53(l are thus additive leading to 

Zf = n Z$ J dR O(R) e^-rW. (61) 
Since this is already first order in n, we only need Eq. l|54Jl up to the zeroth order in n, thus obtaining 

(62) 



where 



Z^( 1 + >„/dRn(R> e «»>) 
pu(R) = pZe J dr'«(r', R)po(r') + ±(/3Ze ) 2 / dr 3 (i> 2 (r, R), (63) 

and hence m extenso 

Z*=Z$(l + \n J dR fi(R) e^ J dr'«(r'3)po(r')+i(/32e ) 2 / dr S (i> 2 (r,R) \ ^ 

Note that in this formula we have the original fugacity A and not the rescaled value A (see Eq. ©). The divergent 
self-interaction is cancelled by an equal and opposite term in /3u(R). The total free energy can now be obtained from 

Z" - 1 

T\ = —k^T lim 



n— >o n 



= § J J drdrVo(r)«(r,r')/Oo(r') + ±Tr 5 (i>(r,r') - fc B T A ^dR fi(R) e^^. (65) 

If we now go back from the grand-canonical to the canonical partition function |24| again via a Legendre transform 
as given in Eqs. I|37[) an d i|38fl . we finally remain with 

T ^ = \ j J drdr'p (r)u(r,r')/9o(r') + ^Tr g{v)v(r,r') - Nk B T\njdR n(R) e' 3 "^. (66) 

The three terms in the above free energy have the following interpretation: the first one is the direct Coulomb 
interaction between the fixed charged surfaces (macroions) of mean charge density pq. The second one, very similar 
to an analogous term in the weak coupling limit, represents the direct effect of the disorder and is proportional to 
the mean-square disorder charge density g. The third term embodies the disorder effect on the strong coupling level: 
It depends in an exponential fashion on the mean-square disorder strength g and has a non-trivial dependence on 
Coulomb interaction potential v(r, r') (see Eq. (jH3J). 

VII. STRONG-COUPLING RESULTS FOR TWO CHARGED PLATES 

Assuming again that our system is composed of two planar surfaces, located at z — ±a, and with surface charge 
density — a, we have 

po(r) = —a5(z — a) — aS(z + a), (67) 

where because of the electroneutrality 

2Sa = NZe a (68) 

(with S being the area for each surface) . The strength of the disorder is assumed to be of the form 

g(r) = gS(z-a)+gS(z + a), (69) 

and the geometry function fl (R) is given by Eq. I|43|l . We can now compute all the terms in the expression for the 
free energy, Eq. I|66|l . explicitly. We obtain 

\ J J drdr'po(r)u(r, r')po(r') = -\ ° S ^ + const. (70) 



11 



where the term const, again stands for all the terms that are independent of the separation between the surfaces. 
Furthermore 

(3u(TL) = [3Ze J dr'v(r', R)po(r') + \{fiZe^f J dr' g(r')v 2 (r', R) = 

a{2a) 1 ((3Zeo) 2 2ng 2 2 
= "/?Ze — ^—t- -T5— ln(a -z ) + const. (71) 

The last term is again independent of the separation between the surfaces and contributes an additive constant to 
the free energy. Inserting all these expressions back into Eq. I|6(jf> . we obtain 

j3T N 1 f3Ze a i m n t( \ 

^T = 2-2^-^ + (x-l)lnU-/(x) 

= i- + (x-l)ln£>-7( X ) (72) 
M 

where D = 2a is the separation between the surfaces, /i = e^j (2-kI-qZg) is just the Gouy-Chapman length for mean 
surface charge density, and 

Z 2 (/3e ) 2 g ^ 7 2„2 n 
8ir(ee y 

is the dimensionless disorder coupling parameter. Note that it is very similar to the electrostatic coupling parameter 
for mean charge 2, Eq. QJ, except that it scales with counterion valency as Z 2 instead of Z 3 . Furthermore 

1 r +D/2 /(D/2) 2 - z 2 \ ~ x/2 
/W = ln ^y D/2 { D 2 ) dZ = 2F ^^ ( ?4 ) 

which is independent of the inter-plate distance D. The effect of disorder on the interactions between disordered 
charged surfaces in the strong-coupling limit is thus contained in the second term in Eq. I|72|l . Its behavior depends 
crucially on whether the dimensionless disorder coupling strength is smaller, equal or larger than one. Evaluating the 
interaction pressure, P(D), from the free energy, Eq. I|72|) . we find, in dimensionless representation, 



where the first term is the standard strong-coupling pressure [24| 

P SC (D) = -1 + ^, (76) 

while the second additive term represents the effects of the disorder (see Fig. EJ. 

From Eq. I|72|l . we can derive also the equilibrium distance D* between the two surfaces, corresponding to zero 
interaction pressure. It is obtained as 

D* = 2 (1 - X ) M- (77) 

In the undisordered case, \ = 0, this reduces to the known form as obtained in Ref. |2^| . that is an equilibrium distance 
equal to twice the Gouy-Chapman length. From this value the equilibrium distance then continuously approaches 
towards zero for increasing \\ it reaches zero for the critical value \ c — 1 and remains at zero thereafter. This 
behavior has all the features of a second-order, quenched-disorder-driven collapse transition with an unusual value of 
the critical exponent (see Fig. [21 inset). 

For x — 1 t ne interaction pressure between the surfaces is obviously constant in the whole range of separations D 
right down to zero. For \ > 1, the counterion entropy, i.e. the — lnD term in Eq. (|72|l . is completely wiped out by 
the surface charge disorder. It thus appears that the quenched charge disorder and the counterion entropy for large 
values of the disorder coupling strength in some sense counteract one another. 
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FIG. 2: Free energy per particle, Eq. I|72|l . as a function of dimensionless inter-plate distance D/(i for \ — 0, 1 and 2 (bottom 
line, middle line, top line), (i = eo/ (27ri?B^o") is the Gouy-Chapman length for mean surface charge density. Inset: dimensionless 
equilibrium distance D* / fi as a function of the disorder coupling parameter \- 



VIII. DISCUSSION 



We have assessed the effect of the quenched macroion charge disorder in the case of interactions between two apposed 
surfaces in a one-component Coulomb fluid. We assumed a Gaussian ansatz for the quenched surface charge density 
distribution, with finite mean and a constant mean-square width. Though the final results within these assumptions 
are very simple they are highly non-trivial and difficult to derive. 

In the case of two ideally polarizable surfaces at separation D, we find that in the weak-coupling or the Poisson- 
Boltzmann level, S <C 1 (where 5 = 2irZ 3 £^(j is the electrostatic coupling due to the mean charge a), the quenched 
surface charge disorder has no effect on the interactions. This is strictly true only for ideally polarizable surfaces. 
In the general case, electrostatic images would contribute to the disorder-driven part of the interaction, making it 
non-zero and dependent on the dielectric mismatch between the dielectric material of the surfaces and that of the 
aqueous solution bathing the electrolyte in between. The disorder-driven part of the interaction in this general case 
can be either attractive or repulsive, depending on the sign of the dielectric mismatch |36) . 

On the strong-coupling level 3> 1, the disorder effects are always present even in the case of ideally polarizable 
bounding surfaces. The disorder-driven part of the interaction free energy turns out to be additive and is linear in the 
disorder coupling parameter x = ^ Z 2 l^g , which is analogous to the mean charge coupling parameter S introduced 
in Ref. |24j , and scales with the square counterion valency Z and the mean-square width of the disorder distribution 
g. For increasing disorder coupling parameter we find that the equilibrium distance is gradually displaced towards 
smaller values if compared to the case with no disorder, x — 0. At a critical value \ c = 1, the equilibrium distance 
becomes zero and remains zero for even larger values of the disorder coupling parameter. Thus one could think of 
the disorder effects in terms of a disorder-driven second-order collapse phase transition. We are unsure at this point 
how general this interpretation might be. In order to assess its validity additional cases of disorder modified Coulomb 
interactions should be investigated. 

Since quenched disorder is ubiquitous in the context of charged biological macroions, our analysis should be of fun- 
damental value especially in this context. In understanding the electrostatic interactions between charged membranes, 
polyelectrolytes and cells, sufficient attention should thus be given to disorder effects that obviously contribute in a 
non-trivial manner to the electrostatic interactions in aqueous electrolyte solutions. Pronounced attraction stemming 
from the strong-coupling electrostatics can thus be made even stronger due to charge disorder effects. 

From our general formulation it follows that even for net-neutral macroions, i.e. a — 0, but with non-zero charge 
disorder, g ^ 0, the total interaction can show electrostatic effects! This is indeed the most non-intuitive and striking 
feature of our theory. Persistent observations of electrophoretic potentials [37j of net uncharged particles might thus 
have an elegant explanation in terms of the disorder effects invoked in this contribution. These might enhance the 
effects of other charging mechanisms such as Lifshitz-van der Waals-interactions driven adsorption or might in fact, if 
proper source of the quenched disorder is identified, be the sole cause of electrostatic effects with nominally uncharged 
macroions. Concurrent measurements of the electrophoretic potentials with an estimation of the local macroion charge 
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disorder could substantiate or disprove the explanation tentatively proposed above. 

Numerical simulations of counterions in the presence of disordered surface charge distributions will be valuable to 
examine the present predictions, in particular, for large mean electrostatic coupling S » 1 and large disorder strength, 
where the disorder-driven collapse transition is predicted. Simulations with disorder effects are however difficult and 
time consuming and we leave them for a subsequent publication [36|. We have also not given appropriate attention 
to the image effects in this contribution. They would add additional details to the disorder-driven interaction that 
would tend to obscure our main points. For reasons of clarity we will thus relegate the image effects to a subsequent 
publication [3(|. 
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